On the apsidal motion of BP Vulpeculae 
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Abstract 



BP Vulpeculae is a bright eclipsing binary system showing apsidal motion. It 
was found in an earlier study that it shows retrograde apsidal motion which 
contradicts theory. In this paper we present the first BV light curve of the 
system and its light curve solution as well as seven new times of the minima 
from the years 1959-1963. This way we could expanded the baseline of the 
O ' investigation to five decades. Based on this longer baseline we concluded that 

^ , the apsidal motion is prograde agreeing with the theoretical expectations and 

^ ' its period is about 365 years and the determined internal structure constant is 

close to the theoretically expected one. 
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CSl ' 1- Introduction 

OO 



The eclipsing n ature of the 10*'' magnitude star BP Vulpeculae was discov- 
ered bv lilies- Almar (1960) an d sh e gave a period value of 1.938 days which was 
. slightly corrected bv lHuthi (Il965l). T hen the star was neglected until the end 



of the 20th century w hen Lacv ( 1992 !) published UBV colours of the system at 
certain phases. Later iLacv et al.l (|2003 ) presented more than 5000 data points 
' in one colour (V) and they have solved their light curve and determined very ac- 

. curately absolute dimensions of the system This F light curve was again solved 

by a si mple, but automa tized code (jPevor . [iooi) and the result was similar to 



that of iLacv et al] (|2003l ). 



The system has an age of about 1 Gyr consisting of A 7 m V -I- F2m V spec 



tral type components and is slightly eccentric (e = 0.0345. iLacv et al.l . l2003f ). 
In eccentric binary star systems the apsidal line is revolving which is called ap- 
sidal motion and therefore the time difference between a primary minimum and 
the subsequent secondary minimum is variable. The apsidal motion is usually 
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caused by two effects: the tidal forces of the components and the effects of gen- 
eral relativity. In the case of a more or less ordinary eccentric binary, i.e. no 
very close third companion in the system, and no extreme stellar rotation, both 
phenomena predict that the apside shows a prograde motion. The total apsidal 
motion is a sum of the classical and relativistic contributions. It is worthy to 
mention that there are several systems, like DI Herculis or AS Camelopardalis 
where the ob served apsida l motions highly differ from the theoretically predicted 
ones (see e.g. lClaretlll998l . and references therein). It seems that the mentioned 



other effects, like third bodies etc. cause these pecularities (for an overview see 



Borkovits et al.L 120071 ) 



Regarding the case of BP Vul, iLacvl (|2003[ ) established a very short apsidal 
motion period (77± 22yr) and they found that it is a retrograde one. If BP Vul 
really would have had a retrograde apsidal motion it would be a new representa- 
tive of the systems which confronts theory and requires further study because a 
retrograde apsidal moti on is i n contr adiction with theory. But the observational 
window in the work of iLacvl (j2003t ) was about one decade only therefore this 
rough estimation for the period of apsidal motion should be refined. For this 
refinement we present BV observations of BP Vul which were obtained 45-49 
years ago and were not published until now. 



2. Observations 

The BV data points of BP Vul published here were obtained by one of us 
(E. lUes-Almar) during the years 1959-1963. The observations were carried out 
by the 60 cm Newtonian-telescope of the Konkoly Observatory, which is located 
at Budapest and it was installed in 1926. The detector was an 1P21 RCA 
photoelectric tube and B and V filters were used. The observations were reduced 
via a standard way. These differential magnitude data (they had been measured 
to a comparison star and transformed into standard system) are published here 
(see Tables 1-2 which are available electronically via the SIMBAD homepag^. 
The B and V curves can be seen in Figure [H 

Using the method of Kwee fc van WoerdenI ( 1956f l from these observations 
we determined seven new times of minima which are reported in Table 3. All 
the minima times we used for the analysis can be found in that Table, too. 



3. Period analysis 

The O — C values were calculated with the following ephemeris: 
Min I = 2 436 860.3311 + 1.9403494 x E 



(1) 



where the initial epoc h was our fi r st pri mary minimum observation while the 
period was taken from Lacv et al. ( 2003h . The O — C diagram is presented in 



^http: / / simbad.u-strasbg.fr/simbad/ 



2 



Table 3: Observed times of minima of BP Vul. Meaning of weights (W): 10: 
CCD/photoelectric minimum, 2: plate minimum, 1: visual observation, 0: not used for the 
calculation because it is an outlier. 



Time of Min. 
(HJD-2 400 000) 


W 


Type 


Ref. 


Time of Min. 
(HJD-2 400 000) 


W 


Type 


Ref. 


21787.723 


2 


P 


1 


37898.427 


2 


P 


2 


22144.750 


2 


P 


1 


37933.365 


2 


P 


2 


22637.569 


2 


P 


1 


37935.294 


2 


P 


2 


23607.741 


2 


P 


1 


38001.256 


2 


P 


2 


25831.362 


2 


P 


2 


38255.437 


2 


P 


2 


26465.843 


2 


P 


1 


38288.425 


10 


P 


3 


26512.439 


2 


P 


2 


38323.326 


2 


P 


2 


26535.718 


2 


P 


1 


38614.433 


2 


P 


2 


26545.469 


2 


P 


2 


38938.446 


1 


P 


4 


26647.308 


2 


s 


2 


38938.452 


1 


P 


4 


26648.332 


2 


P 


2 


38938.453 


1 


P 


4 


26868.522 


2 


s 


2 


41082.509 


1 


P 


5 


26930.399 





p 


2 


41107.549 





P 


5 


27965.784 


2 


p 


1 


41115.516 


1 


P 


5 


28074.410 


2 


p 


2 


41965.387 


1 


P 


6 


28078.299 


2 


p 


2 


42386.456 


2 


P 


1 


28460.527 


2 


p 


1 


42666.772 


2 


s 


1 


29073.744 


2 


p 


1 


43317.780 


2 


P 


1 


29114.476 


2 


p 


2 


43423.529 


2 


s 


1 


29857.637 


2 


p 


1 


44757.540 


1 


P 


7 


31318.694 


2 


p 


1 


44875.875 


2 


P 


1 


33854.692 


2 


p 


1 


45114.551 


1 


P 


8 


34209.792 


2 


p 


1 


45504.566 


1 


P 


9 


34221.446 


2 


p 


2 


45541.452 


1 


P 


10 


34580.427 


2 


p 


2 


45611.295 


1 


P 


11 


35224.585 


2 


p 


2 


45611.302 


1 


P 


11 


35226.576 


2 


p 


2 


45636.490 


2 


P 


1 


35721.355 


2 


p 


2 


45785.882 


2 


P 


1 


36433.446 


2 


p 


2 


45855.778 


2 


P 


1 


36790.450 


2 


p 


2 


45931.461 


1 


p 


12 


36859.3445 


10 


s 


3 


45933.404 


1 


p 


13 


36860.331 


10 


p 


2 


46003.248 


1 


p 


14 


36860.3311 


10 


p 


3 


46290.424 


1 


p 


15 


37116.458 


2 


p 


2 


46321.464 


1 


p 


15 


37438.547 


2 


p 


2 


46356.387 


1 


p 


15 


37506.4677 


10 


p 


3 


46385.508 


2 


p 


1 


37543.3344 


10 


p 


3 


46534.864 


2 


p 


1 


37572.4389 


10 


p 


3 


46612.497 


1 


p 


16 


37642.260 


2 


p 


2 


46612.497 


1 


p 


16 


37867.3714 


10 


p 


3 


46612.500 


1 


p 


16 



1: [Torres fc Guilbaultl (|2003l ): 2: iHuthI (|l965l ): 3: Present paper; 4: BAV 7; 5: BBSAG 
30; 6: BBSAG 12; 7: BAV 34; 8: BBSAG 61; BBSAG 67 10: BAV 38; 11: BBSAG 
69; 12: BAAVSS 61; 13: BBSAG 73; 14: BBSAG 74; 15: BBSAG 78; 16: BRNO 28 
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Table 3: (Continue.) 



Time of Min. 
(HJD-2 400 000) 


W 


Type 


Ref. 


Time of Min. 
(HJD-2 400 000) 


W 


Type 


Ref. 




I 


P 


16 


4Q21 6 450 


I 


P 


29 


4fifil2 502 


1 




16 


4Q21 6 461 


X 


P 


29 


4-fifil 2 50S 


I 


P 


16 


4Q21 6 464 


I 


P 


30 


4-f;fil 2 504 


I 


P 


16 


4Q21 6 468 


I 


P 


30 


46612 506 


1 




16 


4Q21S 393 


I 


P 


29 


46612 507 


1 




16 


49218 402 


1 


P 


29 


46612 50Q 


I 


P 


16 


4Q91 8 41 1 


X 


P 


30 


46614 441 


1 


P 


16 


4921 8 41 3 


I 


P 


31 


46614 445 


1 


P 


16 


49251 390 


1 


P 


30 


4661 4 447 


I 


P 


16 


49321 242 


I 


P 


30 


46614 447 


1 


P 


16 


49934 390 


I 


P 


32 


46614 448 


1 


P 


16 


49967 380 


1 


P 


32 


46614 452 


I 


P 


16 


50002 31 3 


I 


P 


32 


46678 476 


I 


P 


16 


50324 400 


1 


P 


33 


46678 478 


1 


P 


16 


50357 372 


1 


P 


33 


4670Q 547 


2 


P 


I 


50547 551 


I 


P 


34 


46Q7'^ 428 


1 


P 


17 


50681 41 8 


1 


P 


34 


47026 606 


2 


g 


1 


50718 283 


1 


P 


35 


470'HQ S70 


I 


P 


18 


50751 277 


I 


P 


35 


47064 618 


2 


P 


1 


51 036 496 


I 


P 


36 


47'H61 407 


1 


P 


19 


51063 671 7 


10 


P 


37 


47'H6'H 441 


I 


P 


19 


51 1 28 645 


10 




37 


47'^Q2 5S0 


1 


P 


20 


51 1 29 646 


10 


P 


37 


47'^Q2 5S5 


1 


P 


20 


51 327 564 


1 


P 


38 


47392.543 


1 


P 


20 


51364.416 


1 


p 


39 


47396.424 


1 


P 


19 


51397.4114 


10 


P 


40 


47431.331 


1 


P 


21 


51464.3104 


10 


S 


40 


47466.271 


1 


P 


21 


52031.90450 


10 


P 


41 


47788.3674 


10 


P 


22 


52064.89086 


10 


P 


41 


47790.313 


1 


P 


23 


52096.86757 


10 


s 


41 


47823.313 


1 


P 


24 


52098.80834 


10 


s 


41 


48112.405 


1 


P 


25 


52099.8166 


10 


p 


41 


48112.412 


1 


P 


25 


52101.75702 


10 


p 


41 


48147.327 


1 


P 


25 


52164.7794 


10 


s 


41 


48176.432 


1 


P 


25 


52165.78900 


10 


p 


41 


48533.457 


1 


P 


26 


52425.79570 


10 


p 


42 


48723.609 


1 


P 


27 


52487.88765 


10 


p 


42 


48859.422 


1 


P 


28 


52488.81917 


10 


s 


42 



17: BBSAG 84; 18; BBSAG 86; 19: BBSAG 89; 20: BRNO 30; 21: BBSAG 90; 
22: BAV 56 23: BBSAG 92; 24: BBSAG 93; 25: BBSAG 96; 26: BBSAG 99; 27: 
BBSAG 101; 28: BBSAG 102 29: BRNO 31; 30: BBSAG 105; 31: BBSAG 104; 32: 
BB SAG 110; 3 3 : BB SAG 113; 34: BBSAG 115; 35: BBS AG 116 36: BBSAG 118; 
37: iLacv et al.l (|l999l ); 38: BBSA G 120; 39: BRNO 32; 40 lAgerer etHI ^00^; 41: 
iLacv et al.1 (|2002l ): 42: iLacvl l|2002h 
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Table 3: (Continue.) 



Time of Min. 
(HJD-2 400 000) 


W 


Type 


Ref. 


Time of Min. 
(HJD-2 400 000) 


W 


Type 


Ref. 


52495.64880 


10 


P 


42 


53526.9042 


10 


s 


47 


52562.5517 


10 


s 


42 


53527.91289 


10 


P 


47 


52595.5379 


10 


s 


42 


53898.5192 


10 


P 


48 


52724.589 


1 


P 


43 


53933.4432 


10 


P 


49 


52782.8192 


10 


P 


44 


53987.7740 


10 


P 


50 


52814.7964 


10 


s 


44 


54026.5809 


10 


P 


50 


52817.74512 


10 


p 


44 


54325.3939 


10 


P 


51 


53169.8789 


10 


s 


45 


54388.4173 


10 


s 


51 


53186.4111 


10 


p 


46 











43: iDiethelml l|2003l); 44: iLacvl l|2003h : 45: iLacvl (|2004l); 46: IZeidal (|2004l) : 47: 
iLacvl l|2006l '): 48: iHiibscher et al.l l|2006h : 49: iHiibsched (|2007l '): 50: iLacvl l|2007f) : 51: 
IHiibscher et all l|2008D 



Figure [21 It is clear from Figure [2] that the time lag between the primary and 
secondary minima has changed so the apsidal motion is clearly present. In the 
following calculations CCD and photoelectric times of minima had weights of 
10, plate minima had 2, visual observations had 1. 

Using a second-order approximation in the eccentricity, the times of primary 
and secondary minima will occur at the times given below: 

eP 3 e^P 
Min l = To + EP, - — cosw^ H sm2ujE + ■■■ (2) 

TT 8 TT 

P eP 3 e?P 

Min II = To + EPs -cosuje + -- — - sin2w£; - ... (3) 

2 TT 8 TT 

where Tq are the epoch of a primary minimum. Pa is the anomalistic period, Pg is 
the sidereal period, i.e. Pg ~ -Pa(l —^'1'^''^), e is the eccentricity, lue = ujq + uj'E 
where oj' = 2-kP/U . U is the apsidal motion period. 

Applying these formulae we determined the apsidal motion period with 
the upgraded version of LiteAM software developed by T. Borkovits (see e.g. 
Borkovits et al.ll2002l ). We found from the fitting of the 0-C curve that U/P = 



68700±500 and this means U = 365 years, Tq = 51063.6537±0.0001(i7 JD) and 
= 150° ± 5°. This latter value is in good agreement w ith w = 154.7° ± 3.9° 



found by spectroscopic measurements (jLacv et al.l . |2003[ ) . The apsidal motion 
period yields oj « 1.0°/year. 

Note that the O — C curve is not well-covered yet. There is need for more 
observations to determine an exact value of the apsidal motion period in BP Vul 
- our value given above can be regarded as a first approximation. But more 
interesting than the exact value of this period is that the new value yielded a 
prograde motion of the semi-major axis instead of a retrograde one. 
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4. Light curve solution 



F or the light curve solution we used the Wilson-Devinney Code (|Wilsoril . 
The free parameters were the inclination, the dimensionless surface po- 



tentials, argument of periastron and its time-derivative and the luminosities of 
the components. Limb-darkeking coef ficients were fi xed a nd these fixed values 
were interpolated ones from tables of Ivan Hammd (|l993l) . Gravity darkening 
and reflexion coefficients were also fixed. Mass ratio, surface temperatures of 
the components a nd eccentricity of the orbit were fixed at the values given in 



Lacv et al.l (|2003l ). Then differential correction analysis were carried out and 
the stopping criteria was that the change in the parameters in the final step 
should be lower than its standard deviations. Since BP Vul has a fast apsidal 
motion (see previous Section) we used time as an independent variable during 
the modeling rather than phase. The result of the light curve solution can be 
found in Table H 



Comparing our results to the one of lLacv et al.l (|2003f ) we found a remarkably 
excellent agreement in luminosity ratio, but other elements are slightly different. 
However, the precision of our light curve does not reach the precision of their 
one although we have colour i nformation, too. Moreover, they used the so-called 
EBOP code (|Popper &: Etze]| .IT981) which has a slightly different input physics. 
Since we solve these old light curves for the purpose to determine the argument 
of periastron independently, these slight differences do not destroy the validity 
of our light curve solution. 

Thus we concentrate the position of the periastron hereafter. As one can see 
from Table 4 w = 126° ± 5° at epoch HJD =2 436 860.3311 (our adopted epoch) 
which nearly correspond s 1959 October 18. From their spectroscopic measure- 
ments |Lacx!etaD (|2003h had given uj = 154.7° for epoch HJD 2 451023.254 
which nearly corresponds to 1998 July 28. It is easy to compute that these two 
measurements yield 0.74°/yr apsidal motion or U = 486 ± 57 years. 

This value is more than the 365 years apsidal motion period - determined 
from the O — C analysis - by about 120 years. According to us this 33% 
difference is not because the O ~ C diagram is not well-covered and a new light 
curve solution based on multi-colour observations are needed. 

However, the fitted 6j gives 1.0°/yr which is fully in agreement with the 
results of the O — C analysis. Regarding the uncertainty in the position of w 
determined from these old light curves, one can conclude that very likely the 
light curve solution gives lo « 1.0°/yr. 



5. The internal structure constant fc2 

Our next calculations are base d on Gimene^ 1985 ). From the known eccen- 
tricity, masses and period given in Lacv et aLT i 20031 ) one can calculate that the 
relativistic contribution to the apsidal motion in BP Vulpeculae is 7.529 • 10~^ 
degree/cycle. From the observed U — 365 years see above we found a;obs = 
5.24 • 10^'^ degree/cycle. So the Newtonian term in the apsidal motion is 
ioN = ti'obs — LUrc\ = 4.49 • 10^"* degree/cycle. 
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Table 4: Light curve solution of BP Vulpecul ae. Denotions have their usual meaning. Mode 
of the Wilson-Devinney Code was used (seelwiiioEliiii). ri 2 were derived from Qi 2 and 
q by the WD-code itself. 



Quantity 




This paper 


Lacv et al. f2003) 


i 


adjusted 


86.64 ±0.16 


87.67 


Li/Ltot (B) 


adjusted 


0.719 ±0.02 


- 


Li/Ltot (V) 


adjusted 


0.696 ±0.02 


0.718 




adjusted 


126° ± 5° 


154.7° 


Co 


adjusted 


1.006° ± 0.002° 




fii 


adjusted 


7.04 ±0.09 




^1 


adjusted 


6.89 ±0.07 




HJDQ 


fixed 


2 436 860.3311 


2 451 023.254 


mean ri 


derived 


0.162 ±0.027 


0.1931 


mean r2 


derived 


0.141 ±0.031 


0.1552 


9i 


fixed 


1.0 




92 


fixed 


1.0 




Ai 


fixed 


1.0 




A2 


fixed 


1.0 




Ti 


fixed 


7709K 


7709K 


T2 


fixed 


6823K 


6823K 


a^i,boi 


fixed 


0.538 




2^2, bol 


fixed 


0.467 




a;i,B 


fixed 


0.604 






fixed 


0.621 






fixed 


0.534 


0.50 ±0.03 


X2y 


fixed 


0.507 


0.56 ±0.03 


e 


fixed 


0.0345 


0.0345 


9 


fixed 


0.811 


0.811 
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Using the well-known relationship 



C21 + C22 360° 

we found logA:2,obs = —2.66 ± 0.08. Here C21 and C22 are the functions of 
eccentricity, mass ratio and fractional radii of the components and their precise 
form is given in Gimen ez (1985). 

Using the tables of lClaretl (|2004l) with X=0.70, Z=0.02, t=l Gyr and with 



mixing length parameter a = 1.68 and overshooting parameter aov — 0.2 we 
could calculate logfc2,thco — —2.47. Regarding the uncertainties in the deter- 
mined apsidal motion period which appears in the determination of k2 internal 
structure constant we may conclude that the observed and theoretically ex- 
pected values are close to each other. Also note that we have only one secondary 
minimum from the 1960s which is a key point in similar calculations. 



6. Summary 



BP Vulpeculae is an ec centric eclipsing b inary star showing the so-called 
periastron-precession effect. iLacv et al.l (|2003h concluded that this effect causes 
a retrograde motion of the semi-major axis and it has a period of 77 ± 22 
years based on their about one decade long observational material. However, 
retrograde motion contradicts theory. Their explanation was that a possible 
third body in the system could perturb the orbit yielding the observed peculiar 
periastron precession. Nevertheless no spectroscopic evidence was found by 
them for such a third body. 

BP Vul was o bserved at the Ko nkoly Observatory more than forty years 
before the work of Lacv et al] ( 20031 ) by one of the authors of this study. This 
observational material allowed us to determine the times of six primary and 
one secondary minima. With these early observations the baseline could be 
expanded to approximately fi ve decades which w as enough to refine the apsidal 
motion period determined bv iLacv et al.l ( 2003f l. Our O — C analysis based on 
the extended time-line showed that an unseen, dark third body in the system 
cannot be extracted from the presently available minima observations. All of 
these makes very unlikely the presence of a third body with a mass and orbit 
which would cause a peculiar periastron precession. 

First time two-colour light curves were presented by us for BP Vul. The light 
curve solutioi i - using the Wilson-De yinney Cod e - yie lded very similar results 
comparing to Devor ( 2005[ )'s one and Lacv et al. (2003)'s one . In addition, the 
O — C analysis in this paper showed that the apsidal motion period in BP Vul 
is prograde and it has a period of about 365 years - it is in agreement with 
the value determined with less accuracy from the light curve. The prograde 
motion means that BP Vul is not a representative of problematic cases and it 
is in agreement with theoretical expectations. Ne vertheless we conc luded that 
the negative apsidal motion rate determined by iLacv et aL (2003) is only a 
consequence of their short observational window. 
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Figure 1: The B (top) and V (bottom) light curves of BP Vul obtained in the years 1959-1963. 
Bottom is the B curve while top curve is the V one which is shifted by 0.3 magnitudes for the 
sake of clarity. 

We also calculated the k2 internal structure of BP Vul and found it being 
close to the theoretical value. The slight difference should be refined in the 
future with a better observed O — C diagram without gaps. Therefore the 
minima observations of BP Vulpeculae in the future are needed. 

The comments on the first version of the manuscript by Drs J. Jurcsik and 
K. Olah is acknowledged. 
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Figure 2: The O — C diagram of BP Vulpcculac. Filled and oijcii squares represent primary 
and secondary minima, respectively. Lines show the fits to the O — C residuals for the primary 
and secondary minima, respectively. The weights of the different kind of minima can be found 
in the text. 
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